Electronic structure and magnetic properties of NaOsC>3 
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A comprehensive investigation of the electronic and magnetic properties of NaOsOa has been made 
using the first principle calculations, in order to understand the importance of Coulomb interaction, 
spin-orbit coupling and magnetic order in its temperature-induced and magnetic-related metal- 
insulator transition. It is found that its electronic structure near the Fermi energy is dominated by 
strongly hybridized Os 5d and O 2p states. Despite of the large strength of spin-orbit coupling, it 
has only small effect on the electronic and magnetic properties of NaOs03. On the other hand, the 
on-site Coulomb repulsion affects the band structure significantly, but, a reasonable U alone cannot 
open a band gap. Its magnetism is itinerant, and the magnetic configuration plays an important 
role in determining the electronic structure. Its ground state is of a G-type antiferromagnet, and 
it is the combined effect of U and magnetic configuration that results in the insulating behavior of 
NaOsOa . 

PACS numbers: 71.20.-b, 71.30.+h, 72.80.Ga 
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I. INTRODUCTION 

It is well known that the Coulomb interaction 
among 3c? electrons in transition-metal oxides (TMO) 
is substantially important, which induces peculiar 
properties, such as metal-insulator transition^, colos- 
sal magnetoresistance 2 - and high critical temperature 
superconductivity 3 . Because the 5 c? orbitals are highly 
extended compared to those in the 3c? systems, it is nat- 
ural to expect that the electronic correlations are weak 
and have only negligible effect in the 5 c? compounds. 
However, recent theoretical and experimental works have 
given the evidence on the importance of Coulomb interac- 
tions here^ 5 -. On the other hand, the spin-orbit coupling 
(SOC) in the 5c? transition metal elements is expected to 
be strong due to the large atomic number™. Hence, due 
to the interplay of electron correlations and strong spin- 
orbit interactions, various anomalous electronic proper- 
ties have been observed/proposed in the 5 c? transition ox- 
ides, such as 3 e ff — 1/2 Mott stated, giant magnetoelec- 
tric effect, high T c superconductivity 8 -, Weyl semimetal 
with Fermi arcs 9 -, Axion insulator with large magneto- 
electric coupling-™ topological insulato r 11 ' 12 , correlated 
metali 3 -, Kitaev modeM, etc. 

One class of the well studied 5 c? compounds are the 
osmatesi£~— . For example, the physical properties of 
Cd 2 0s 2 07 are quite intriguing. It has been found that 
Cd20s207 is metallic at room temperatures, while un- 
dergoing a metal-insulator transition (MIT) at about 
230 K 15 . Experiments reveal that this MIT is contin- 
uous and purely electronic. Moreover, it is coincident 
with a magnetic transition of antiferromagnetic (AFM) 
characte r 15 : 16 . Therefore, experimentalists argue that 
Cd20s207 is the first well-documented example of a pure 
Slater transitio n 16 : 26 . However, despite of the vast efforts 
devotecU 5 - - — , its exact magnetic ground state configura- 
tion is still unknown due to the strong geometric frus- 
tration of the pyrochlore lattice. Therefore, theoretical 



evidence of the Slater transition in this compound is still 
lacking. 

Recently, using high pressure technique, Shi et al£L 
synthesized another osmate: NaOs03. Similar to 
Cd20s207, NaOsOa also exhibits a temperature-induced 
MIT, which is again accompanied by a magnetic ordering 
without any lattice distortion 2 ^-. However, better than 
Cd 2 0s 2 07, NaOs0 3 has a simple perovskite structure, 
consequently being free from the complication induced 
by magnetic frustration. Therefore, NaOsOa provides a 
unique platform to understand the temperature-induced 
and magnetic-related MIT. Based on the experimental 
crystal structure, Shi et al^ also perform the band- 
structure calculation for this compound. They 2 ^ find 
that both LDA and LDA+SO calculation give the para- 
magnetic solution. Their numerical results 21 show that 
Coulomb U is not efficient, and antiferromagnetic corre- 
lation is essential to open the band gap2i. Recently, there 
are experimental and theoretical evidences of the impor- 
tance of electronic correlation and spin-orbital coupling 
in 5d transition-metal compounds. Therefore a compre- 
hensive investigation of the effect of the Coulomb interac- 
tion, SOC and magnetic order on its electronic structure 
and MIT is still an interesting problem which we address 
in the present work. 



II. METHOD 

The electronic band structure calculations have been 
carried out by using the full potential linearized aug- 
mented plane wave method as implemented in WIEN2K 
package 27 . Local spin density approximation (LSDA) 
is widely used for various 4c? and 5c? transition metal 
oxides^ ' 25 ' 28 : 29 , and we therefore adopt it as the 
exchange-correlation potential. The muffin-tin radii for 
Na, Os and O are set to 1.13, 1.02, and 0.90 A, respec- 
tively. The basic functions are expanded to R m tK ma2; =7 
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(where R mt is the smallest of the muffin-tin sphere radii 
and K ma2 . is the largest reciprocal lattice vector used in 
the plane wave expansion), corresponding to 1915 LAPW 
functions at the T point. Using the second-order varia- 
tional procedure 3 ^, we include the spin-orbital coupling 
interaction (SOC), which has been found to play an im- 
portant role in the 5d system^^iS. A 10x6x10 mesh is 
used for the Brillouin zone integral. The self-consistent 
calculations are considered to be converged when the dif- 
ference in the total energy of the crystal does not exceed 
0.1 mRy and that in the total electronic charge does not 
exceed 10~ 3 electronic charge at consecutive steps. 

TABLE I: Numerical and experimental internal coordinates 
ofNaOsOa. 





calculation 


experiment 


atom 


x y z 


x y z 


Na 
01 
02 


0.0392 1/4 0.9910 
0.4919 1/4 0.0885 
0.2940 0.0428 0.7046 


0.0328 1/4 -0.0065 
0.4834 1/4 0.0808 
0.2881 0.0394 0.7112 



NaOsOs has an orthorhombic perovskite structure 
with space group of Primal. There are four formula 
units (f.u.) per unit cell, and the 20 atoms in the unit cell 
can be classified as four nonequivalent crystallographic 
sites: Na, Os, Ol and 02 according to the symmetry. 
They are located at 4c, 46, 4c and 8d sites, respectively 
and result in seven internal coordinates. From the X-ray 
diffraction experimental, the lattice constants of NaOsOs 
are determined to be a=5.384 A, &=7.580 A and c=5.328 
A, respectively. Based on the experimental lattice pa- 
rameters, we optimize all independent internal atomic 
coordinates until the corresponding forces are less than 
1 mRy/a.u. We confirm that the Coulomb U and SOC 
have only small effect on the crystal structure and list in 
Table I the internal atomic coordinates by LDA calcu- 
lations. Our numerical internal coordinates are in good 
agreement with the experimental result, as shown in Ta- 
ble 1. 



III. RESULTS AND DISCUSSIONS 

Using the experimental lattice constants and the nu- 
merical internal coordinates, we first perform non- 
magnetic LDA calculation. The total density of states 
(TDOS), Os 5d partial density of states (PDOS), Ol 2p 
PDOS, 02 2p PDOS, Na 2s, 2p PDOS has been plotted 
in Fig.l(a)-(e), respectively. Our TDOS is very similar 
to that obtained based on the experimental crystal struc- 
ture (See Fig. 4a of ReL^i). The energy range, -9.0 to -2.4 
eV is dominated by Ol 2p and 02 2p bands with a small 
contribution from Os 5d state. Both Na 2s and 2p states, 
appearing mainly above 4 eV, have also considerable dis- 
tribution between -9.0 to -2.4 eV, where O 2p state is 
mainly located, indicating the non-negligible hybridiza- 
tion between Na and O states despite that Na is highly 
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FIG. 1: Electronic density of states (DOS) from LDA calcu- 
lation. Fermi energy E/ is set to zero, (a) TDOS, (b) Os 5d 
PDOS, (c) Ol 2p PDOS, (d) 02 2p PDOS, (e) Na PDOS. 




FIG. 2: Band structure of NaOsOs, shown along the high 
symmetry directions, (a) LDA, (b) LDA+SO. 



ionic. The Os atom is octahedrally coordinated by six O 
atoms, making the Os 5d band to split into the t2 g and e g 
states, and the 12 ta g bands are located from -2.8 to 1.2 
eV, as shown in Fig. 2a. Due to the extended nature of 
5d states, the crystal splitting between t2 S and e g states 
is large, and the e g states are located about 2.0 eV higher 
than the Fermi energy (Ey) and disperse widely. While 
providing the basic features of the electronic structure, 
LDA produces a metallic state due to partially occupied 
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FIG. 3: Electronic density of states (DOS) from LDA+SO+U 
(U=2.0 eV) calculation. Fermi energy Ey is set to zero, (a) 
FM configuration, (b) G-AFM (c) C-AFM, (d) A-AFM (e) 
D-AFM. The black/red line is for up/down-spin, respectively. 
For AFM configurations, spin-up is the same as spin-down, 
thus only one spin channel is plotted. 



Os 5d t2 g band. 

SOC of bd electrons is about 0.5 eV— , which is one or- 
der of magnitude larger than that of 3d electrons. There- 
fore, SOC usually changes the 5d band dispersion signif- 
icantly and plays an essential role in the gap opening of 
Si'2lr04 as well as of pyrochlore iridates^^. In order to 
investigate the effect of SOC on the electronic structure, 
we compare the results obtained in the presence and ab- 
sence of SOC, which are given in Fig. 2. The difference 
between the bands with and without SOC is small, as 
demonstrated in Fig. 2. For 5d 5 electronic configuration 
of S^IrO^, and A2L2O7 (^4=Y or rare earth)2, where 
3 e ff — 1/2 picture is valid, SOC has a dramatic effect 
on the band structure. In NaOs03, Os occurs in its 5 + 
valence and there are 3 electrons in its t2 9 band. Since 
t2g band is half filled, it is natural to expect the effect of 
SOC to be small. As shown in Fig. la, the Fermi level 
is located near a sharp peak in the DOS. The relatively 
high density of states at the Fermi energy (N(E/)) sug- 
gests the possibility of a Stoner instability against ferro- 
magnetism (FM) . Therefore, we perform a spin polarized 
calculation, however our LSDA+SO calculation with ini- 
tial FM setup converges to the non-magnetic state. Thus 
consistent with Shi et al31, the FM state is not stable in 
LSDA+SO calculation. 



TABLE II: Spin (S) and orbital (O) moment (in fi B ) as well 
as the total energy E to t per unit cell (in eV) for several mag- 
netic configurations, as calculated using LDA+U+SO method 
with U=2.0 eV. {Etot is defined relative to the G-AFM con- 
figuration.). 

Configuration G-AFM FM C-AFM A-AFM D-AFM 
E tot 0.243 0.186 0.282 0.205 

(5") 0.94 0.22 0.54 0.29 0.20 

(O) -0.11 -0.01 -0.04 -0.03 -0.03 



Although the 5d orbitals are spatially extended, it has 
been found that the electronic correlations are impor- 
tant for 5d TMO^£. Moreover, the experiment reveal 
that NaOsOa has a long-range magnetic order at low 
temperature 2 ^. We therefore utilize LSDA+U scheme^, 
which is adequate for searching for magnetically ordered 
insulating ground states^ 2 -. Although, the accurate value 
of U is not known for perovskite osmates, the estimates 
of the values of U have been recently obtained between 
1.4 and 2.4 eV in layered Sr 2 Ir04/Ba 2 Ir04 33 . We gen- 
erally expect the screening to be larger in three dimen- 
sional (3D) systems than in two dimensional (2D) ones, 
and one can image that U in NaOsOs should be smaller 
than that in Sr 2 Ir04/Ba 2 Ir04. We therefore perform 
LDA+U+SO calculation and vary parameter U between 
0.5 and 2.0 eV. Numerical results show that the electronic 
correlations can stabilize the FM configuration and nar- 
row the Os ta g band. However, as shown in Fig. 3(a), 
our LDA+U+SO calculation with U= 2.0 eV still gives 
a metallic solution. Naively, one may expect that using 
larger Coulomb U will result in an insulating state. How- 
ever, consistent with Shi et al<21, our additional calcula- 
tions show that increasing U cannot solve this problem, 
and even a quite large U (=6.0 eV) cannot open the band 
gap. Therefore, electronic correlations alone cannot ex- 
plain the insulating behavior, and the MIT is not of a 
Mott-type. 

After studying the effect of SOC and U, we subse- 
quently investigate the effect of various magnetic orders. 
We considered four antiferromagnetic (AFM) configu- 
rations besides the FM state: A-type AFM state (A- 
AFM) with layers of Os ions coupled ferromagnetically 
in a given set of (001) planes but with alternate planes 
having opposite spin orientation; C-type AFM state (C- 
AFM) with lines of Os ions coupled ferromagnetically in 
a given direction (001) but with alternate lines having 
opposite spin orientation; G-type AFM state (G-AFM) 
with Os ions coupled antifcrromagnetically with all of 
their nearest neighbors; D-type AFM state (D-AFM) 
where Os ions lying within alternating planes perpen- 
dicular to [001] direction are coupled ferromagnetically 
along either [010] or [100] directions while different lines 
are coupled antiferromagnetically. Same as with the FM 
setup, the LSDA+SO calculation with U—0 for all con- 
sidered AFM setups converges to the nonmagnetic metal- 
lic state. 

On the other hand, the non-zero Coulomb interaction 
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U of Os 5d is found to stabilize the AFM configura- 
tion. Our calculation confirms that the magnetic order 
has a significant effect, and for a reasonable U (<2.0 
eV), G-AFM configuration is the only insulating solu- 
tion as shown in Fig. 3. Moreover, regardless the value of 
U, the G-AFM configuration always has the lowest to- 
tal energy. Thus we believe G-AFM configuration is the 
magnetic ordering state observed by the experimental . 
With increasing U the band structure will change, but 
only when U is larger than 1.0 eV, the G-AFM solution 
becomes insulating. The DOS from U=1.0 eV (see the 
Fig. 4c of Ref.21) is similar with that from U=2.0 eV (see 
Fig. 3b of present work), which again indicates that the 
Coulomb U is not efficient to open the band gap. It is 
found that the magnetic moment is mainly located at 
Os site, and despite of strong hybridization between Os 
5d and O 2p, O site is basically non-magnetic (less than 
0.003 /tig). The numerical data for U =2.0 eV are given in 
Table II. For the 5d 5 electronic systems such as Balr03, 
Sr2ir04, pyrochlore iridates etc, it has been found that 
due to the strong spin-orbit entanglement in 5d states, 
the magnetic orbital moment is about twice larger than 
the spin moment^^, even in the presence of strong 
crystal field and band effects. Contrary to 5d 5 systems, 
the obtained orbital moment for NaOsOa is much smaller 
than its spin moment, showing again that SOC effect 
is small for this 5d 3 electronic configuration case. As 
shown in Table II, the magnitude of magnetic moment 
is sensitive to the magnetic configuration, indicating the 
itinerant nature of magnetism. For the same U value, 
the G-AFM configuration always has the largest mag- 
netic moment among the considered states. However, as 
shown in Table II, our numerical magnetic moment (0.83 
/j, B ) is much smaller than the experimental one^i. For 
an itinerant magnet, one may still fit the x(^) curve by 
the Curie-Weiss law, but cannot estimate the magnetic 
moment accurately based on the Curie- Weiss constant 35 . 
Thus, the experimental magnetic moment may not be re- 
liable. The energy difference between various magnetic 
configurations is large, which is consistent with the ob- 
served high magnetic transition temperature (about 410 



K)2i although here one cannot estimate the interatomic 
exchange interaction and T n based on the difference be- 
tween total energies accurately as in the local moment 
systems 3 ^. Since the G-AFM configuration is the only 
insulating state, it is easy to understand that both mag- 
netic and electronic phase transitions occur at the same 
temperature and our calculation indeed confirms that the 
MIT of NaOsOa is a Slater-type transition. 



IV. SUMMARY 

In summary, we have investigated the detailed elec- 
tronic structure and magnetic properties of NaOsOs 
using full potential linearized augmented plane wave 
method. Our results show that the electronic structure 
near the Fermi energy Ej is dominated by strongly hy- 
bridized Os 5c? and O 2p states. Despite of its big value 
the SOC has only weak effect on the band structure and 
magnetic moment. The electronic correlations alone can- 
not open the band gap, and the low temperature phase 
of NaOs03 is not a Mott-type insulator. The magnetic 
configuration has an important effect on the conductiv- 
ity, and the ground state is a G-type AFM insulator. 
It is the interplay of the Coulomb interaction and mag- 
netic ordering that result in the insulating behavior of 
NaOs0 3 . 
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